#ifdef _TEST_F_
#include "Spline.h"

void complete_spline1()
{
    INTERPC interpC;
    interpC.points = {1, 2, 5, 6, 7, 8, 10, 13, 17};
    interpC.interpSites[1] = {3.0, 1.0};
    interpC.interpSites[2] = {3.7};
    interpC.interpSites[5] = {3.9};
    interpC.interpSites[6] = {4.2};
    interpC.interpSites[7] = {5.7};
    interpC.interpSites[8] = {6.6};
    interpC.interpSites[10] = {7.1};
    interpC.interpSites[13] = {6.7};
    interpC.interpSites[17] = {4.5, -0.67};

    Spline<1, 3, ppForm> spline = interpolate<3, ppForm>(interpC, complete);
    spline.plot('b');
}

void complete_spline2()
{
    INTERPC interpC;
    interpC.points = {17, 20, 23, 24, 25, 27, 27.7};
    interpC.interpSites[17] = {4.5, 3.0};
    interpC.interpSites[20] = {7.0};
    interpC.interpSites[23] = {6.1};
    interpC.interpSites[24] = {5.6};
    interpC.interpSites[25] = {5.8};
    interpC.interpSites[27] = {5.4};
    interpC.interpSites[27.7] = {4.1, -4.0};

    Spline<1, 3, ppForm> spline = interpolate<3, ppForm>(interpC, complete);
    spline.plot('b');
}

void complete_spline3()
{
    INTERPC interpC;
    interpC.points = {27.7, 28, 29, 30};
    interpC.interpSites[27.7] = {4.1, 0.33};
    interpC.interpSites[28] = {4.3};
    interpC.interpSites[29] = {4.1};
    interpC.interpSites[30] = {3.0, -1.5};

    Spline<1, 3, ppForm> spline = interpolate<3, ppForm>(interpC, complete);
    spline.plot('b');
}

int main()
{
    complete_spline1();
    complete_spline2();
    complete_spline3();

    return 0;
}

#endif
